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ABSTRACT 

Soft- inplane rotor systems can significantly reduce 
the inplane rotor loads generated during the maneu- 
vers of large tiltrotors, thereby reducing the strength 
requirements and the associated structural weight of 
the hub. Soft-inplane rotor systems, however, are 
subject to instabilities associated with ground reso- 
nance, and for tiltrotors this instability has increased 
complexity as compared to a conventional helicopter. 
Researchers at Langley Research Center and Bell 
Helicopter- Textron, Inc. have completed an initial 
study of a soft-inplane gimballed tiltrotor model sub- 
ject to ground resonance conditions in hover. Para- 
metric variations of the rotor collective pitch and 
blade root damping, and their associated effects on 
the model stability were examined. Also considered 
in the study was the effectiveness of an active swash- 
plate and a generalized predictive control (GPC) al- 
gorithm for stability augmentation of the ground res- 
onance conditions. Results of this study show that 
the ground resonance behavior of a gimballed soft- 
inplane tiltrotor can be significantly different from 
that of a classical soft-inplane helicopter rotor. The 
GPO-based active swashplate was successfully imple- 
mented, and served to significant ly augment damping 
of the critical modes to an acceptable value. 


This paper is declared a work of the U.S. Government and 
is not subject to copyright protection in the United States. 
Presented at the 42nd AIA A/ASME/ASCE/AHS/ASO 
Sturctures, Structural Dynamics, and Materials Conference, 
Seattle, Washington. April 16-19, 2001. 


INTRODUCTION 

While three-bladed, gimballed. stiff-inplane rotor 
systems are employed for the current generation of 
tiltrotors (XV-15, V-22, BA609), the weight and per- 
formance penalties of stiff- inplane rotors may become 
too significant for use on future systems larger than 
the V-22. In airplane mode, rotor loads are usually 
small compared to those in helicopter mode, due to 
reduced disk loading and axisymmetric inflow con- 
ditions. However, during airplane-mode maneuvers 
sucii as a rapid pull-up, high oscillatory inplane ro- 
tor loads are developed. Tiltrotors employ large 
highly-twisted blades, creating a significant aerody- 
namic forcing in the inplane direction where the cen- 
trifugal restoring forces are much lower than those 
associated with the flapping direction. StifF-inplane 
hubs must be designed for these inplane loads us- 
ing additional structure which increases weight and 
can make the inplane stiffness higher still, leading to 
more increases in the inplane loads. For larger tiltro- 
tors, this cyclic design challenge can lead to very high 
structural weights or even to an infeasible design al- 
together. As a design trade-off, the maximum aero- 
dynamic load capability of a large tiltrotor may be 
limited to reduce the structural weight of the hub. 
Such is the case for the V-22 which currently em- 
ploys a controller to limit the body pitch-rate motion 
of the aircraft and thereby curtail rotor inplane loads 
which could otherwise rise above their design allow- 
ables (Ref. 1). 

Soft-inplane rotor systems can significantly reduce 
the inplane rotor loads generated during maneuvers 
of larger tiltrotors. thereby reducing strength re- 
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Figur : Effect of lag frequency on hub loads. 

quirements on the hub, leading to reduced structural 
weight and improved aircraft agility. The plot of 
Fig. 1 illustrates the benefit with respect to loads 
for moving from a stiff- inplane design point to a soft- 
inplane design point. For a soft-inplane rotor system, 
the lower the inplane natural frequency becomes, the 
lower the associated hub loads become. The plot 
also shows that there is a theoretical asymptotic in- 
crease in the load factors at IP due to the resonance 
of the inplane mode (lead-lag mode) with rotor speed 
( “lag mode” will be used throughout the remainder of 
the paper as a description of this mode, however, the 
terms “inplane” and “lag” will be used interchange- 
ably to be consistent with common usage of the in- 
dustry). For a soft-inplane rotor system to be ef- 
fective, it should have a lag frequency below about 
0.9/rev, otherwise the loads are about the same as 
those of the stiff inplane rotor system. A desirable 
design point for development of a soft-inplane tiltro- 
tor is with a lag frequency of about 0.7/rev (as in- 
dicated on Fig. 1). This lag frequency provides a 
good reduction in design loads and it is about the 
lowest lag frequency that is practical for the types of 
soft-inplane hubs currently under consideration for 
full-scale development. Soft- inplane rotor systems, 
however, are subject to instabilities associated with 
ground resonance. Before soft-inplane rotor systems 
become viable for application to tiltrotors, a com- 
prehensive understanding of the means for avoiding 
ground resonance instabilities must be established, 
but there have been few investigations related to this 
subject to date. 

The ground resonance phenomenon is well- 
understood for conventional helicopters, and the re- 


search presented in Refs. 2-10 explain well the phe- 
nomenon and characteristics of the onset of ground 
resonance instabilities. However, there are several 
differences between helicopters and tiltrotors which 
suggest that ground resonance behavior of tiltrotors 
may be significantly different than ground resonance 
of a helicopter. There is increased complexity for 
a tiltrotor associated with the participation of elastic 
wing modes, in addition to the body and rotor modes, 
which couple through the changing rotor speed during 
wind-up (start-up). The ground resonance condition 
is improved with an increase in damping in both the 
rotor system lag mode and the body or wing mode 
with which it participates. While the damping of the 
body modes may be supplemented through mechan- 
ical design of the landing gear, there is little that 
can be adjusted in the elastic modes of a tiltrotor 
wing without difficulty and significant wing design 
changes. 

The type of hub used for a tiltrotor is another 
important consideration influencing the characteris- 
tics of ground resonance. Previous studies, as is 
discussed in the review paper of Ref. 10, have con- 
sidered the important differences between the fully- 
articulated, hingeless, and bearingless rotor systems 
with respect to ground resonance on helicopter sys- 
tems. However, these past efforts may not be ap- 
plicable to the present effort because none of these 
conventional hub systems are currently employed on 
tiltrotors. Standard bearingless rotor systems do not 
provide enough pitch control for tiltrotor applications 
while hingeless rotors do not provide enough flap] ; 
control during differential cyclic maneuvers. Arti<. i- 
lated rotor systems are heavier and produce higher 
drag than the gimballed rotor systems that have be- 
come the standard for tiltrotors to date. The current 
study will focus on the ground resonance behavior 
of a soft- inplane gimballed rotor system, w ich may 
display results significantly different from those asso- 
ciated with the hub systems considered in past stud- 
ies. However, it should be noted that this type of 
rotor system may not be feasible for implementation 
on tiltrotors due to other stability considerations such 
as high-speed whirlflutter. Although fully-articulated 
rotors may produce a greater drag penalty in airplane 
mode, some modification of this hub type may need to 
be considered to make soft-inplane tiltrotors feasible 
for all aeroelastic considerations. Additional analy- 
sis and testing is required to identify the ideal hub 
configuration for soft-inplane tiltrotor applications. 

To date, the most significant studies of tiltro- 
tor systems subject to ground resonance conditions 
were made by Boeing Helicopter with the develop- 
ment of their Model 222 tiltrotor concept, a hinge- 
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less soft-inplane rotor system, which was their en- 
try into a NAS A/ Army-sponsored tiltrotor research 
aircraft program (eventually named the XV-15 and 
won by the Bell Helicopter stiff-inplane rotor sys- 
tem). Ground and air resonance behavior of the 
Boeing soft- inplane configuration was addressed in 
several experimental and analytical studies using dif- 
ferent size rotor test apparatuses, beginning with a 
1/ 10-scale wind-tunnel model as described in Ref. 11, 
and ending with a full-scale 26- ft. diameter semispan 
model tested in the NASA Ames 40- x 80- ft. tun- 
nel as described in Refs. 12 and 13. The Boeing 
soft-inplane design had a relatively high inplane nat- 
ural frequency, such that the design rotor speed in 
hover mode did not create a ground resonance prob- 
lem. The only experimental results associated with 
aeromechanical instability of this configuration was 
obtained with the system in airplane mode subject 
to air resonance conditions. This configuration was 
also not gimballed, it was a standard hingeless soft- 
inplane design. Thus, the tiltrotor configuration 
of the current study is unique in two important de- 
sign parameters: 1) the use of a gimbal and constant 
velocity joint which has a significant effect on the 
Coriolis forces and therefore the flap-lag blade cou- 
pling, and 2) the use of a “low” lag frequency (about 
0.5/rev) which creates a resonance condition between 
the low-frequency lag mode and the critical tiltrotor 
pylon/ wing mode that is well within the design ro- 
tor speed envelope. While this latter parameter is 
not considered a. desirable design goal, it does pro- 
vide some benefits for the current study as shall be 
discussed. 

Researchers at NASA Langley Research Center 
(LaRC) and Bell Helicopter- Textron, Inc. (BHTI) 
have completed an initial study of a soft-inplane gim- 
balled tiltrotor model subject to ground resonance 
conditions in hover. This effort was planned as part 
of a Memorandum of Agreement between the NASA 
LaRC and BHTI to “Perform Experimental Aeroelas- 
tic Studies of a Tiltrotor Model,” and represents the 
first experiment in a series of tests to be conducted 
with this soft-inplane hubs on the Wing and Rotor 
Aeroelastic Testing System (WRATS). The phrase 
“initial study”, as used in the first sentence of this 
paragraph, indicates that the current hub design is 
not scaled from a full-scale design that is deemed 
practical for large scale application. Rather, the cur- 
rent design was developed as a low-cost soft-inplane 
modification to the existing stiff-inplane gimballed 
hub representative of the V-22. The objective of 
the current study was to evaluate the stability char- 
acteristics of the gimballed soft-inplane rotor system 
subject to ground resonance conditions in hover, and 



Figure 2: Hover-cell home of the Wing and Rotor 
Aeroelastic Testing System at the TDT. 

also to assess the feasibility of using a generalized pre- 
dictive controller (GPC) with an active swashplate 
control system to augment damping and eliminate 
the ground resonance instabilities. Parametric vari- 
ations of the rotor collective pitch and blade root 
damping, and their associated effects on the model 
stability were also examined in the study. 

APPARATUS 

The experimental study was performed in a 30’ x 
30 1 hover cell (Fig. 2) located in a high-bay building 
adjacent to the Transonic Dynamics Tunnel (TDT) 
at the NASA Langley Research Center in Hamp- 
ton. Virginia. Notable features of the facility are 
a backstop mount which centers the tiltrotor model 
in the hover cell at the same height above the floor as 
when mounted in the TDT wind-tunnel test section 
(8 ft.); a snubber stand to halt pylon motion in the 
advent of an instability; a 3000 psi hydraulic sys- 
tem; 100 psi shop air; 440 volt electric motor power; 
and a closed-circuit chilled-water system for motor 
cooling. The model mount was designed to provide 
stiffness similar to that of the TDT test section side- 
wall mount so that system frequencies are identical 
in either the hover or wind-tunnel facilities. Instru- 
mentation wiring runs from the model into an un- 
derground cable tray, out of the hover cage, and into 
a block house control room. Testing is monitored in 
the control room using a closed-circuit camera system 
and standard televisions. Signal conditioning, data 
acquisition equipment and the model pilot control 
console are located side-bv-sidc in the control room. 
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i igure 3: WRATS girnballed soft-inplane parametri- 
cally variable hub. 

A remote control unit is used to operate the MG set 
which controls model rotor speed. It should be noted 
that while several comparisons between the hover cell 
and TDT test section were mentioned in this para- 
graph, the current Study w r as performed entirely in 
the hover test cell. No wind-tunnel test results have 
yet been conducted on the soft-inplane version of the 
til trot or model. 

The hover test facility described in the previous 
paragraph was developed in 1995 as a dedicated area 
for the Wing and Rotor Aeroelastic Testing System 
(WRATS), a semi-span 1 /5-size aeroelastic tilt rotor 
model based on the V-22. This tiltrotor model has 
been used in several aeroelastic experimental efforts 
beginning in 1984 as part of the Navy’s JVX pro- 
gram, and more recently has been on loan to NASA 
LaRC. Since 1994, BHT1 and NASA LaRC have had 
an ongoing cooperative research agreement in place 
to perform experimental aeroelastic investigations in- 
volving the WRATS model with several associated 
modifications. Some of the more notable investi- 
gations include stability augmentation using a com- 
posite tailored wing (Ref. 14), vibration reduction 
using an active flap (Ref. 15). and vibration reduc- 
tion using an active swashplate (Ref. 16). Important 
general features of the model are listed as follows: 
an aeroelastically-scaled wing with removable airfoil 
panels, a dynamically-scaled pylon with a downstop 
spring tuned to provide elastic mode shapes and fre- 
quencies close to those associated with the full-scale 
conversion actuator (different springs are used for dif- 
ferent conversion actuator positions), a girnballed 3- 
bladed hub with a constant-velocity joint, and a set 
of aeroelastically-scaled rotor blades. The current 
study was performed using the WRATS model, but 


Figure 4: WRATS soft-inplane system in hover. 

with the normal stiff-inplane hub modified with the 
addition of lag hinges to provide a soft-inplane ro- 
tor configuration. As illustrated in Fig. 3, the soft- 
inplane hub is parametrically variable with a set of re- 
placeable coil springs (located behind the lag spring 
catches) used to tune the lag stiffness, and an ad- 
justable hydraulic damper used to tune the lag damp- 
ing. The dampers are adjusted by turning a screw 
which controls the flow of fluid within the damper 
casing. With the screw open to about 5 turns the 
damper is ineffective, and only the structural and 
mechanical damping of the system is present (about 
3%). With the screw shut, the damper is almost 
immovable and the lag motion is critically damped. 

Another notable feature of the WRATS tiltro- 
tor model is the hydraulically-controlled swashplate 
which has high bandwidth controller capability. 
Three oil cylinder actuators are used to position the 
swashplate, each controlled by a Moog servo valve 
using an attached linear variable displacement trans- 
ducer (LVDT) for position feedback. The pilot con- 
trol console has three inputs for AC/DC signals (ac- 
tive control commands) which may be summed with 
the pilot DC input commands and sent to the swash- 
plate control system. The model is shown in Fig. 4 
configured for hover testing; and surrounded by the 
snubber system, but loosely attached to it so that the 
system properties are not effected. 

For testing ground resonance instability, the snub- 
ber system provided two benefits. In addition to pro- 
viding a means for preventing damage to the model 
with the onset of an instability, it also provided 
a mechanism for performing isolated rotor testing. 
While the snubber was activated to restrain the up- 
per portion of the pylon, the lower portion of the 
pylon could be clamped to the support stand. This 




Figure 5: WRATS tiltrotor model setup for isolated 
rotor testing. 

combination would raise the lowest fixed-system fre- 
quencies to about 20 Hz. The isolated rotor test 
setup was used to determine frequencies and damp- 
ing for the lag mode in the uncoupled system, which 
are beneficial to know' before initiating the part of the 
experiment w'here instabilities are likely to occur. A 
photo of the model setup for isolated rotor testing, 
with snubber activated and lower pylon clamped to 
the support stand, is shown in Fig. 5 . 

PARAMETRIC STUDY IN HOVER 

The stability characteristics of a gimballed soft- 
inplane tiltrotor model were determined experimen- 
tally as a function of collective pitch and lag damper 
settings for a lag frequency setting of 0.5/rev. The 
choice of the lag frequency used for the study is 
a crucial one because it defines the rotor speed at 
which ground resonance will occur. In general, the 
lower the lag frequency, the greater the problem with 
ground resonance instability. While a practical 
tiltrotor configuration may most likely have a lag fre- 
quency of 0.7/rev to 0.8/rev, the current test was de- 
veloped using a much lower value as indicated. The 
reason for using 0.5/rev was two-fold: 1) there is a 
limited number of stock coil springs that could be 
modified to fit into the spring cage which would pro- 
duce reasonable lag frequencies, and 2) this value 
of lag frequency provided instabilities at low rotor 
speeds so that the ground resonance related charac- 
teristics of the system could be fully examined be- 
fore, during, and beyond the rotor speed at which 
resonance occurred. 

The system frequencies and associated damping 
were determined experimentally by exciting the py- 


lon with a cyclic stick-stir at the natural frequency of 
the mode of interest. After removing the excitation, 
data w r as processed for approximately 5 seconds, and 
the frequency and damping were determined using 
a moving block method (Ref. 17) from the digitized 
time histories. Obtaining the lag response required 
slightly more work at some rotor speeds because the 
response was highly damped and thus did not pro- 
vide enough cycles outside of the “noise” level to 
accurately calculate a frequency or damping of the 
mode. To determine the frequency of the lead lag 
mode in these instances, the amplitude of the lag re- 
sponse was measured during the excitation produced 
by the cyclic stick stir. This measurement was re- 
peated for several perturbations of the excitation fre- 
quency. When the excitation matched up to the 
natural frequency of the lag mode then the response 
amplitude of this mode reached a peak. 

The fixed-system mode of interest for the tiltro- 
tor model in hover is best described as a wing tor- 
sion/chord coupled elastic mode at approximately 4.9 
Hz. For hover, the pylon is oriented vertically above 
the wing elastic axis creating a significant mass-offset 
which creates the large coupling between the chord 
and torsion modes. A second wing mode, beam (ver- 
tical) bending, occurs at 5.1 Hz, but since the beam 
motion is not highly-coupled with the lag rotor mo- 
tion (the dominant associated pylori motion is per- 
pendicular to the rotor plane), there is a negligible 
participation of this mode in the ground resonance 
phenomenon of interest. Also, it should be noted 
that while in airplane mode the wing beam mode is 
highly coupled with the wing torsion mode, this is 
not the case for the hover mode, because of the po- 
sition of the pylon and the location of its associated 
mass offset. A third mode at 11.5 Hz is best de- 
scribed as a wing chord /torsion mode, but because 
of the relatively high frequency, this mode is also in- 
significant with respect to the ground resonance con- 
dition. The 4.9 Hz wing torsion/chord mode (here- 
after referred to as WTC) and rotor lag frequencies 
are shown as a function of rotor speed in Fig. 6 for a 
condition of 8° collective pitch and damper setting of 
3% critical (nonrotating). The plot shows that there 
is only a slight drop in the wing frequency with rotor 
speed. The regressive low'- frequency lag mode be- 
comes a progressive mode at about 490 RPM. As the 
rotor speed continues to increase the low-frequency 
lag mode and the WTC mode become more coupled 
and the associated frequencies begin to converge. 

The damping of these modes is illustrated in the 
plot of Fig. 7. The plot show's that the lag mode is 
highly damped, greater than 3% critical except over 
the rotor speed range w'here an instability is develop- 
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Figure 6: Coupled frequencies of the wing and rotor 
modes of interest for ground resonance (8° collective 
and 3% critical lag damper setting). 

ing. Damping at very low rotor speeds could not be 
measured because there is not sufficient aerodynamic 
forcing to produce excitations using the stick stirs. 
The WTC mode damping is shown to decrease at a 
nearly constant rate with rotor speed, and the mode 
eventually becomes unstable at about 720 RPM. 

The variation of WTC mode damping with rotor 
speed is shown in Fig. 8 for three collective pitch set- 
tings. The results show that the system behavior is 
extremely sensitive to collective pitch, and this trend 
is more sensitive than that associated with a typi- 
cal helicopter system. With an increase in collective 
pitch, the aerodynamic forces at higher rotor speeds 
may produce significant hub forces which couple the 
WTC and blade lag modes, leading to an aeroelastic 
instability. The WTC mode damping shown in Fig. 8 
as a function of rotor speed at the 5° and 12° collec- 
tive pitch settings also illustrates a trend which is un- 
like that of classical helicopter ground resonance, but 
more like a helicopter pylon instability. In classical 
ground resonance conditions with insufficient rotor 
lag damping, a sudden and steep reduction followed 
by a recovery in the lag mode damping is often ob- 
served in the vicinity of the resonant frequency. This 
plot shows that the 5° and 12° collective cases have 
a relatively slow, but constant, reduction in WTC 
fixed-system damping as rotor speed increases. Also, 
following the instability, the WTC mode damping 
does not recover at higher rotor speeds. 

The parametric effect of the damper setting was 
also considered in this study. Surprisingly, there 
was little difference in the results of the fixed-system 



Figure 7: Damping of the wing and rotor modes of 
interest for ground resonance (8° collective and 3% 
critical lag damper setting). 

damping and the associated ground resonance insta- 
bility, even with significant changes in the lag damper 
setting. In fact, these results are not shown using a 
plot, as all the curves essentially are overlaid with 
those shown in Fig. 8. While there is suspicion that 
these results may be a consequence of the particular 
design, and that some structural nonlinearities as- 
sociated with the lag mechanics may play a role in 
this behavior, this remains a significant and surpris- 
ing result of the study. It can also be seen that the 
damper does have an effect on the lag motion itself. 
The plot of Fig. 9 shows the IP response of the lag 
motion as a function of rotor speed, in the vicinity 
of the lag mode IP cross-over. The amplitude of 
the response is shown to be greatly influenced by the 
damper setting, but this effect does not seem to in- 
fluence stability of the system. 


GENERALIZED PREDICTIVE CONTROL 
Control Algorithm 

This section describes the theory of GPC in gen- 
eral terms. Implementation of GPC to the specific 
model used in the study is described in a subsequent 
section. The essential features of the GPC adaptive 
control process are depicted in Fig. 10. The system 
(plant) has r number of control inputs u, m number 
of measured outputs t/, and is subject to unknown 
external disturbances d. Measurement noise is also 
present. There are two fundamental steps involved: 
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Rotor Speed, RPM 

Figure 8: Effect of collective pitch on ground reso- 
nance stability with lag damper set to 3% critical. 

(1) identification of the system; and (2) use of the 
identified model to design a controller. A finite- 
difference model in the form of an auto- regressive 
moving average with exogenous input (ARX) model 
is used here. This model is used for both system 
identification (SID) and controller design. System 
identification is done on-line in the presence of any 
disturbances acting on the system, as indicated in the 
center box of the diagram. In this way, an estimate 
of the disturbance model is reflected in the identi- 
fied system model, and does not have to be modeled 
separately. This approach represents a case of feed- 
back with embedded feedforward. Because the dis- 
turbance information is embedded in the feedforward 
control parameter, there is no need for measurement 
of the disturbance signal (Ref. 18). The parameters 
of the identified model are used to compute the pre- 
dictive control law. A random excitation u i( j (some- 
times called dither) is applied initially with the vector 
of closed-loop control inputs u c equal to zero to iden- 
tify the open-loop system. Dither is added to u c if it 
is necessary to re-identify the system while operating 
in the closed-loop mode. 


The relationship between the input and output 
time histories of a MIMO system are described by the 
time-domain auto-regressive exogenous (ARX) finite- 
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Figure 9: Effect of lag damper on IP lead-lag re- 
sponse (8° collective). 

difference model: 

V(k) = 

<*iy(k - 1) -f ot< 2 y(k — 2) + ... + a p y(k - p) 

+ fto u (k) + fi\u(k — 1) + ... + l3 p u(k - p) (1) 

This equation states that the current output y(k) at 
time step k may be estimated by using p sets of the 
previous output and input measurements, y(k - 1), 
..., y(k — p) and u(k — 1), .... u(k — p); and the cur- 
rent input measurement u(k). The integer p is called 
the order of the ARX model. The coefficient matri- 
ces a 7 and /?,; appearing in Eqn. 1 are referred to as 
observer Markov parameters (OMP), or ARX para- 
meters, and are the quantities to be determined by 
the identification algorithm. 

Closed-loop robustness is enhanced bv performing 
the system identification in the presence of the exter- 
nal disturbances acting on the system, thereby ensur- 
ing that disturbance information will be incorporated 
into the system model. The goal of SID is to deter- 
mine the OMP based on input and output data. The 
OMP may be determined by any SID techniques that 
return an ARX model of the system. 

The ARX model is used to design the controller and 
leads to a control law that in the case of a regulator 
problem has the general form given by: 

ttr(fc) = 

a c ]y (k - 1) + a$y{k - 2) + ... + a«y(k - p) 

4* >3 ^ u(k) -f- (i^iik — 1) + ... + if p u(k — p) (2) 



Rotor Speed, RPM 




This equation indicates that the current control in- 
put u c (k) may be computed using p sets of the previ- 
ous input and output measurements. The coefficient 
matrices and 0^ appearing in Eqn. 2 are the con- 
trol gain matrices. 

System identification in the presence of the oper- 
ation disturbances acting on the system is the first 
of the two major computational steps. The external 
disturbances acting on the system are assumed to be 
unknown (unmeasurable). The number of control in- 
puts is r and the number of measured outputs is ra. 
The system is excited with band-limited white noise 
for SID. These random excitations are input to all r 
control inputs simultaneously and the corresponding 
rn responses are measured. The input and output 
time histories ( u and y) are digitized according to 
the l number of time points, and these parameters 
are then used to form the data matrices y and V as 

V = YV (3) 

where 


y = [y(0) y(l) 2 /( 2 ) ... yip) - y(l - 1)] (0 

and 



u(0) u(l) u(2) . 

• u(p) 

.. u(l- 1 ) 


t»(0) w(l) • 

■ v(p- 1 ) 

... v(l- 2) 

v = 

?'(<•) • 

• v(p - 2 ) 

... v(l — 3) 



»(0) 

. . . v(l — p — 1 ) 


(5) 

The order of the ARX model (p) and the number of 
time steps (/) must be specified by the user. The size 
of y is m x l and the size of V is [r + (r + m)p] x /. 
Equations 4 and 5 follow from writing the discrete- 
time state-space equations for a linear time-invariant 
system at a sequence of time steps k = 0, l, .... (/ — 1) 
and grouping them into matrix form. The vector v 
is defined as 

''<*>={$)} ,6 > 

which has size (r + m) x 1. In forming the matrices 
of Eqns. 4 and 5, it has been assumed that the state 
matrix A is asymptotically stable so that for some 
sufficiently large p, A k ~ 0 for all time steps k > p, 
and that an observer has been added to the system. 
It is through these expedients that the matrix V is 
reduced to a size amenable for practical numerical 
computation of its pseudo- inverse. The SID process 
yields OMP rather than system Markov parameters 
(SMP) because of the inclusion of an observer. A 



Figure 10: The GPC adaptive control process block 
diagram. 

complete discussion of these aspects of the develop- 
ment may be found in Ref. 19. 

Y is the matrix of observer Markov parameters 
which are to be identified and has the form 

Y = { 0 Q 0 1 ai 0 2 a 2 0 3 a 3 ... 0 V a,, } (7) 

By solving Eqn. 3 for Y, the solution may be written 
as 

Y = yV~ x = yV T [VV T )~ l ( 8 ) 

If the product VV T is well-conditioned, the ordinary 
inverse may be taken. Otherwise, a pseudo-inverse 
must be used. It should be noted that because the 
size of VV T is much smaller then V , a pseudo-inverse 
may be appropriate even if the product is well condi- 
tioned. 

At this point in the derivation the output predic- 
tion equation has been established for one time step. 
The control law is enhanced by developing the output 
prediction equation several time steps ahead (multi- 
step approach). With the a* and determined by 
the system identification process, the multi-step out- 
put prediction equation may be derived at time step 
k -b j. Equation 1 may be written at time step k + j 
as 


y(k + j) = 

a[ j) y(k - 11 •• a 2 j) y(k - 2) + ... + a ( p j) y(k - p) 

+ 13q u(k + j) + 0q u{k + j — 1) + ••• + 0lpu(k) 

+ 0[ j) u(k - 1) + 0)pu{k - 2) + ... + 0pu(k - p) 

(9) 

This equation shows that the output y(k ; + j) at time 
step k + j may be estimated by using p sets of the 
previous output and input measurements, y(k- 1), 
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y(k — p) and u(k — 1), ..., u(k — p), and the (unknown) 
current and future inputs, u(fc), u(k+ 1), ..., u(k+ j). 
The user selects a discrete value for the prediction 
horizon h p which provides the range of values for j 
to be considered. Letting j in Eqn. 9 range over 
the set of values j = 1, 2, .... h p — 1, the resulting 
equations can be assembled into a multi-step output 
prediction equation: 

Vh p (k) = Tu he (k) + Bu p (k - p) + Ay p (k - p) (10) 

where the coefficient matrices T, B, A are formed 
from combinations of the observer Markov parame- 
ters cti and f3 i . The quantity yh p {k) is the vec- 
tor containing the predicted future outputs, whereas 
u hc (k) is the vector containing the future control in- 
put yet to be determined. The quantities u p (k — p) 
and y p (k - p) are vectors containing the previous p 
sets of control inputs and outputs, respectively. The 
vector Uh c in the derivation of Eqn. 10 is originally 
the size of h p , but by introducing a new user se- 
lectable quantity defined as the control horizon h Cj 
where h r < h p , the controller performance may be 
improved. The objective of the controller is now to 
predict the output for h p time steps ahead, given the 
input for h c steps ahead. 

The goal of the GPC algorithm is to determine the 
set of future commands u(k), u(k + 1), ..., u(k+h c -l) 
that are required to achieve a desired predicted re- 
sponse y(k), y(k + 1), y(k + h p — 1). The multi- 
step output prediction equation, Eqn. 10, is used to 
define an objective function whose minimization with 
respect to Uh t . (k) leads to the control law from which 
a vector of future control inputs can be computed us- 
ing the p sets of previous control inputs and measured 
outputs. The predictive control law is obtained by 
minimizing the deviation of the predicted controlled 
response (as computed from the multi-step output 
prediction equation) from a specified target response 
over a prediction horizon h p . To this end, an error 
function s is defined which is the difference between 
the desired (target) response yr(k) and the predicted 
response y hp ( k): 

£ = yrik) - yh t> (k) = 

yr(k) - Tu hc (k) - Bu p (k - p) - Ay p (k - p) (11) 

An objective function J is then defined which is 
quadratic in the error e and in the unknown future 
controls : 

J = e T Re + uJ lc Qu} lc (12) 

As shown, two weighting matrices are introduced: Q 
is used to limit the control effort and stabilize the 


closed-loop system, and R is used to weight the rela- 
tive importance of the differences between the target 
and predicted responses. Both matrices are sym- 
metric and positive- semidefinite, and typically are as- 
sumed to be diagonal. Also, it is typical in practice 
for Q to have the same value w c along its diagonal 
and for R to have the same value w r along its diago- 
nal. Minimizing J with respect to Uh ( (A') and solving 
for the same yields: 

Uh c { k) = 

- (T t RT+ Q)~ l x 

T T R{-y T (k) + Bu p (k ~ p) -f Ay p (k - p)] (13) 

which is the control sequence to be applied to the 
system over the next h r time steps. The first r val- 
ues (corresponding to the first future time step) are 
applied to the r control inputs. The remainder are 
discarded, and a new control sequence is calculated 
at the next time step: 

Ur(k) = 

- ? f yr(k) -f (f'u p (k - p) + c\ ( y p (k - p) ( 14) 

The target response is zero for a regulator problem 
and non-zero for a tracking problem. Matrix Q must 
be tuned to ensure a stable closed-loop system, and 
typically h c is chosen equal to h p . However, h c may 
be chosen less than h p resulting in a more stable, but 
sluggish, regulator. 

GPC Implementation 

For the ground resonance investigation, active con- 
trol of the WRATS tilt rotor model was performed 
using the high-bandwidth capability (50 Hz) of the 
swashplate control system (3 hydraulic actuators). 
The GPC commands were summed at the pilot con- 
trol console with the pilot input commands before be- 
ing sent to each of the three actuator control valves. 
The GPC computations were distributed between two 
CPUs on a single computer with no user interaction 
required for data transfer. On user command, data 
required for system identification was collected on 
CPU1 and sent to CPLT2 where SID was performed 
and the control gain matrices, a <: and (?\ were com- 
puted. The control parameters were automatically 
sent to CPU1 which used the p latest data sets to 
compute the control commands (?q.) to be sent to 
the sw r ashplate actuators. If re- identification was re- 
quired, the process was repeated on user command. 
All the algorithms were implemented in dSPACE (a 
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Figure 11: Effect of the GPC controller on ground 
resonance related damping. 

commercially available real-time control system) on 
a PC with 500 MHz CPUs. The adaptive sequence 
described in the GPC algorithm section was not in- 
cluded in the current implementation of GPC. 

Key results of the investigation are illustrated in 
the plot of Fig. 11. For the baseline system (GPC 
off) and a collective pitch setting of 8°, the WTC 
mode damping is shown to decrease until the model 
becomes unstable at slightly over 700 RPM. With 
the GPC control system activated, a large increase 
in damping was immediately apparent. The plot 
shows a consistent increase of about 2% critical, and 
the system did not become unstable within the rotor 
design speeds considered. Although no results are 
shown, the vibratory loads associated with the wing 
modes were examined with GPC both on and off. A 
reduction in vibratory loads of about 50% was a gen- 
eral rule for the settings on the controller that were 
optimized for the stability augmentation on the sys- 
tem. It also should be noted that a great deal of work 
was performed in obtaining the settings (/, p, /i r , h p , 
w c , and w r ) that produced both controller stability 
and significant damping increases for the model. It 
was not uncommon during the process of tuning the 
controller to develop a control instability and then 
have to shut the control system down and/or snub 
the model. Once the control settings were estab- 
lished, the performance of the GPC controller were 
impressive as indicated by Fig. 11. 


CONCLUDING REMARKS 

With respect to ground resonance behavior, the ex- 
perimental results obtained for the current test of 
a gimballed-hub tiltrotor configuration were signifi- 
cantly different than those generally associated with 
classical helicopter soft-inplane rotor systems. In 
classical ground resonance conditions with insuffi- 
cient rotor lag damping, a sudden and steep reduction 
followed by a recovery in the lag mode damping is of- 
ten observed in the vicinity of the resonant frequency. 
The results of the current study show a relatively 
slow, but constant, reduction in the WTC (critical 
body mode) damping as rotor speed increases. Also, 
following the instability, the WTC mode damping 
does not recover at higher rotor speeds. As there 
does appear to be an aeroelastic nature to the insta- 
bilities in hover, the sensitivity of the instabilities ob- 
served to aerodynamic-related design variables such 
as pitch-flap coupling and pitch-lag coupling will need 
to be studied further. 

Testing of a new generalized predictive control 
(GPC) system was also completed in this study. The 
use of GPC with an active swashplate proved to be 
highly effective at increasing damping and eliminat- 
ing ground resonance induced instabilities. In ad- 
dition to increasing the damping, the active control 
system also reduced vibratory response by about 50% 
simultaneously. While these results are encouraging, 
there was some difficulty encountered in determining 
controller settings that were both effective and stable. 
The active control implementation used was also not 
completely automated nor adaptive. Current efforts 
are focused on improving these aspects of the GPC 
controller. 
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